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Abstract. A new method for measuring the shear induced by 
gravitational light deflection is proposed. It is based on ana- 
lyzing the anisotropy induced in the auto-correlation function 
of the extragalactic background light which is produced by 
very faint distant galaxies. The auto-correlation function can 
be measured 'locally', and its anisotropy is caused by the tidal 
gravitational field of the deflecting mass distribution in the 
foreground of these faint background galaxies. 

Since the method does not require individual galaxy detec- 
tion, it can be used to measure the shear of extremely faint 
galaxies which are not detectable individually, but are present 
in the noise. The shear estimated from the auto-correlation 
function of the noise provides an independent measurement 
which can be compared to the shear obtained from the dis- 
tortion of individual galaxy images. Combining these two in- 
dependent estimates clearly increases the sensitivity of shear 
measurements. In addition, our new method may allow to de- 
termine the local magnification caused by the deflector if the 
auto-correlation function is caused by a large number density of 
faint galaxies; in this case, the intrinsic auto-correlation func- 
tion may provide a 'standard source' with respect to which 
shear and magnification can be obtained. Applications to real 
and synthetic data are shown and the feasibility of our new 
method is demonstrated. In particular, we present the shear 
maps obtained with our new method for the double QSO 
2345+007 and the cluster C10024+16 and compare them to 
published shear maps. 
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1. Introduction 

Weak gravitational lensing is one of the most promising tools 
for probing the amount and the distribution of mass in the uni- 
verse. The small anisotropy in the shape of galaxies induced 
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by gravitational shear maps almost directly the total mass pro- 
jected near the line of sight to these faint galaxies regardless 
of its dynamical state. Recent theoretical and observational 
work on galaxy-galaxy lensing (Brainerd et al., 1996, Schnei- 
der & Rix 1996), weak lensing in rich lensing-clusters (Kaiser & 
Squires 1993, Seitz & Schneider 1995, Bonnet et al. 1994, here- 
after BMF; Fahlmann et al. 1994, Squires et al. 1996, Seitz et 
al. 1996), and from large scale structures (Blandford et al. 1991, 
Kaiser 1992, Villumsen 1996) justifies the expectation that out- 
standing results on the mass distributions of these systems will 
be available in a near future which will solve numerous puz- 
zling issues on the (dark) mass distribution, on the number 
density of dark halos, and on the relation between light and 
mass of gravitational structures. 

In the work quoted above, the shear is inferred from the el- 
lipticity measured on individual background galaxies (we shall 
refer to this method in the following as 'standard method', 
to contrast it to our new method). In the case of weak lens- 
ing, the induced distortion is extremely small which makes the 
measurement of the ellipticity from ground-based telescopes 
images a major concern. The measurement and the various 
intrumental and atmospheric corrections of shape parameters 
have been discussed by Kaiser, Squires & Broadhurst (1995) 
and Bonnet & Mellier (1995, hereafter, BM). From these inves- 
tigations, one finds that dramatic limitations come from rather 
poor signal-to-noise ratios of images at very weak shear level. 
Ultra-deep CCD images on wide fields allow the computation 
of the averaged ellipticity on a large area with many galaxies 
which increases notably the signal-to-noise ratio. However, this 
strategy has important shortcomings: 

1. The CCD area over which ellipticities can be averaged 
to obtain a local estimate of the shear is limited to a typical 
scale across which the magnification matrix is approximately 
constant. Hence, it depends on the typical angular scale of 
the gravitational structure. Furthermore, averaging on large 
areas decreases the spatial resolution of the shear map, which 
eventually washes out shear produced by substructures. 

2. This limitation can be overcome by ultra-deep images 
which provide a large number of galaxies on small areas. Un- 
fortunately, because these objects are so faint, their centroid 
and their shape parameters can be determined only poorly 
from observation. Furthermore, deep HST images show that 
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a large fraction of the faintest galaxies have no well-defined 
shapes and show many bright spots which can be mistakingly 
considered as galaxy centroids on ground based images. 

In this paper we propose an alternative method for the 
shear measurement which can overcome these critical issues. 
Such a method, which uses the auto-correlation function of 
the brightness distribution (hereafter ACF) of the whole im- 
age, was implemented by some authors (Cole et al. 1992, Dal- 
canton 1996) to detect small amplitude variations in the mean 
noise associated with faint cluster of galaxies. The philosophy 
of the correlation methods is to avoid individual object de- 
tection in order to work at the noise limit, and consequently 
to go much deeper than classical analyses. Here, the funda- 
mental point is that a stretching of a shear-free image by 
gravitational distortion induces an equivalent stretching in the 
auto-correlation function. The strong advantage is that the co- 
addition is made without computation of the central moments 
of individual objects, and, as shown later, the signal-to-noise 
of the ACF can be greater than a direct co-addition of the ob- 
jects. In fact, the concept of object itself is no more relevant: 
the CCD field is viewed as a fluctuating 'density field', and the 
objects are viewed as 'overdensities'. Consequently, any over- 
density present on the CCD, even if it is related to extremely 
faint sources which have not been detected with standard de- 
tection algorithms, provides information about the ACF. For 
that reason, the faintest sources, which only emerge from the 
noise limit of the background, will be detected by the ACF. It 
will provide additional information of the shear field from the 
modification of the noise pattern which, as a result, appears as 
a "sheared noise pattern" . 

Obviously, the signal-to-noise ratio of the ACF depends 
on the number density of faint background sources. The most 
recent deep galaxy counts in the visible at the Keck tele- 
scope (Smail et al. 1995) yields compelling evidence that the 
ACF should be considerably more powerful than the standard 
method. Consider for instance a standard detection criterium 
based on the positive measurement of a connex-shape object 
formed with at least 6 pixels at least one a (standard devia- 
tion) above the background. In practice, this definition rejects 
smaller objects with flux per pixel larger than two a above 
the background because they look like cosmic ray events (steep 
slope of the contours). Therefore, the largest "undetected" pat- 
tern which could be a galaxy in the CCD may have 5 pix- 
els, having at most one a above the background. From our 
own observational data obtained at the Canada- France-Hawaii 
Telescope (CFHT), this "undetected" pattern corresponds to 
objects with magnitude V = 30.4 on a 5 hours exposure. Ex- 
trapolating the deep galaxy counts of Smail et al. (1995) we 
expect about 2 xlO gal./deg /mag. at V = 30.4, or about 550 
gal./arcmin 2 /mag.! In our CCD field this amounts to about 
one object per field of 5 arcsec 2 ! The ACF method will take 
into account any of these ultra-faint galaxies without the need 
to determine their individual position and shape. We therefore 
expect a considerable increase of the signal. Conversely, we 
may use the increase of the signal-to-noise ratio on the shear 
to analyse the projected number density of these ultra-faint 
sources. 

In this first paper we only discuss the implication of the 
ACF technique for the (weak) shear measurement. The analy- 
sis of the background population is discussed in a forthcoming 
paper (Van Waerbeke et al.; Paper II). In the next section, 
we present the basic formalism and demonstrate how the ACF 



is "gravitationally distorted" by the shear. We show that the 
ACF method is strictly equivalent to a single object analysis 
but with a high signal-to-noise ratio. The third section shows 
how to extract the gravitational shear information from the 
ACF. In particular, we present a detailed analysis of the error 
estimation. The fourth section validates the method on simu- 
lated data and on real images from which we know the shear 
pattern from the standard detection technique (BMF). There, 
we compare shear maps obtained from the auto-correlation 
function with those measured by the 'standard method' for 
the double QSO 2345+007 and for the inner part and an outer 
region of the cluster CI 0024+16 and show that these shear 
maps are in excellent agreement. 

2. The Auto-correlation Method 

2.1. The basic formalism 

Throughout this paper we keep solely a geometrical point of 
view of the shear effects. This means that we do not consider 
the mass reconstruction problem. 

Consider a gravitational lens at some redshift Zd which is 
sufficiently small such that the distance of all sources can be 
approximately set to a constant (or more precisely, the distance 
ratio D dB /D B is the same for all sources). In a solid angle which 
is sufficiently small such that the local lensing properties do 
not change over this region, the lens equation can be locally 
linearized, such that the angle (3 of a light ray in the 'source 
plane' is related to its observed angular position in the lens 
plane by 

(3 = f3 + A(0 )(0-e ) (1) 

where Oo is a central position of the angular region considered, 
/3 () the corresponding angular position in the source plane, and 
A is the local Jacobian matrix, which reads 

M 1- -;; 71 i-7+ 71 H^-^) <*> 

where k is the local dimensionless surface mass density, 7,, 
i = l,2, are the two components of the shear (for more details, 
see Schneider, Ehlers & Falco 1992), X is the two-dimensional 
identity matrix, and 

, . / cos 2^5 sin 2ip \ , . 

J ^=Un2^ -cos 2 J ® 

is a symmetric traceless tensor which describes the orientation 
ip of the shear. Conservation of surface brightness implies that 
the observed surface brightness 1(0) is related to the intrinsic 
one 7 (s) (/3) by 

I{6) = &\A0) (4) 

where we have applied a translation of the coordinates in the 
source plane to remove the constant terms in (1); this can be 
done since absolute positions in the source plane are unobserv- 
able anyway. 

We now define the auto-correlation function of the bright- 
ness in the image plane by 

£(0) = <(/(#)-/)(/(# + 0)-7)> t? (5) 
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where the angular brackets denote an average over all positions 
and 7 is the mean surface brightness of the image. Inserting 
the relation (4) into (5), one finds immediately that 



C(0)=C (s) (^) 



(6) 



where the auto-correlation function ^ is defined in analogy 
to (5) with I replaced by I^ a \ Hence we see that the 'observed' 
correlation function £ is related in a simple way to the 'intrin- 
sic' one. In particular, the intrinsic correlation function can 
be assumed to be isotropic; then (6) shows that the observed 
correlation function will in general be anisotropic. The power 
spectrum of the surface brightness, defined as the square of the 
Fourier-transform £(k) of the correlation function, is related to 
the intrinsic power spectrum by 



|ak)| 2 =M 2 |£ (s) 04- 1 k)| 



(7) 



where p = 1/detA is the local magnification of the lens. 

If the intrinsic correlation function £^(/3) is known, the 
relation (6) allows to determine the local Jacobian matrix A, 
and thus the shear components 7 and surface mass density k 
(note that this is not strictly true if the lensing is not 'weak', 
as will be discussed later). We shall return to the practical 
questions concerning this issue further below. However, even if 
the intrinsic correlation function is unknown, but only assumed 
to be isotropic (which should be a very good approximation), 
one can obtain some information on the Jacobian matrix A, as 
we show next. 

Consider the tensor of second moments of the auto- 
correlation function, 
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J d 2 6>g(6>) OjOj 



(8) 



f si 

and the corresponding definition for the tensor M\- of second 
moments of the intrinsic auto-correlation function; using the 
relation (6), one readily obtains 



Ma = A^ 1 AJ1 1 M$ 



(9) 



an equation very similar to the transformation of second bright- 
ness moments of individual galaxy images (e.g., Blandford et 
al. 1991, Schneider & Seitz 1995, BM). The main formal dif- 
ference is that individual galaxies cannot be assumed to be 
isotropic in the source plane, so that the observed ellipticity of 
a galaxy image depends on both, the Jacobian matrix and the 

(s) 

intrinsic source ellipticity. In the case considered here, m:/ 
can be considered an isotropic tensor, so that 



MSi 



(10) 



where we have defined the mean quadratic radius M of the 
intrinsic correlation function in the second step. Inserting (10) 
into (9) and using the form (2) of the matrix A, one finds after 
some algebra 



M = 



M 



(1 



2\2 



-*) 2 (i-M 2 ) 
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V 2 ff2 1 
where we have defined the reduced shear 



(11) 



9i ■= 



7i 



(12) 



(1-kY 

We shall see that it is convenient to also introduce the distor- 
tion 

2g t 



5, := 



(1 + lfll 2 )' 



(13) 



which has also been used in weak lensing studies before (e.g., 
Miralda-Escude 1991, Schneider & Seitz 1995, BM). With this 
definition, (11) can also be written as 



M = 



M(i + \ g \ 2 ) n+di 82 \ 
(l-^(l-\g\ 2 ) 2 ^ l - 5 ^ 



(14) 



From the latter form, one can immediately read off that 



So = 



Mu - Mi 

tr~M 
IMii 



(15) 



tiM ' 

where XxM = Mu + M22 is the trace of M. Furthermore, 
from the trace of (14), one obtains 



(I--) 2 



M 



4det7M 



(trM±2VdetM). 



(16) 



In order to appreciate this solution, the corresponding situ- 
ation for the determination of lensing parameters from indi- 
vidual galaxy images should be recalled: there, using only the 
shape of the galaxy images, together with the assumption that 
the intrinsic orientation of galaxies is randomly distributed, 
the distortion S can be uniquely determined. Furthermore, the 
local magnification can in principle be obtained, either from 
changes of the local number density of background galaxy im- 
ages due to the magnification (anti-)bias (Broadhurst, Taylor 
& Peacock 1995), or by changes of the image sizes at fixed 
surface brightness (Bartelmann & Narayan 1995). The same 
information can be obtained in the situation under considera- 
tion here: the distortion can be uniquely determined (15) from 
assuming that the intrinsic correlation function is isotropic, 
without knowledge of its radial dependence, whereas the (ab- 
solute value of the) magnification 



VdetM 
M 



(17) 



can be determined if the second moment of the intrinsic corre- 
lation M function is known. Note that a given combination of 
<5 and \n\ does not specify the shear 7 and local surface mass 
density k uniquely; in fact, for each such combination, four 
different pairs of (7, k) are obained, corresponding to the four 
different combinations of (/i > 0, fi < 0) and (k < 1,/t > 1). 
Only one of these combinations, fj, > 0, k < 1, is relevant for 
non-critical clusters, or the region outside the critical curves in 
critical clusters. 

We can now consider the case of weak lensing, i.e., k«1, 
7 <IC 1. Assuming again that the intrinsic correlation function 
is isotropic, £ (s) (/3) =£ (s) (l/3|), we can expand (6) to obtain 



e(0)«£«(0)-0£W'(0) [K + 7 Jfo,-0)], 



(18) 



where we have written — #(cos 1?, sin $) . Equation (18) shows, 
in the case of weak lensing, that the observed correlation func- 
tion is the sum of the intrinsic correlation function and a per- 
turbative term which linearly depends on both k and 7. 
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2.2. Inclusion of a point spread function 

Equation (4) is only an idealized description of the brightness 
distribution obtained from observations; in practice, it is af- 
fected by observational effects (seeing, tracking errors, etc.) 
and noise. We consider the first of these effects here, and as- 
sume that it can be described by a point spread function (PSF) 
W (ip) which is taken to be normalized, 

J d 2 i>W(tp) = 1. (19) 
The observed brightness profile /( obs )(<?) is then given by 
/(°bs)(0) = J d 2 <pW(9-ll>)I(lp) 

= ( d> W{6 - ip) I (B) {Alp). (20) 



This observed brightness profile can now be inserted into 
the definition (5) to obtain the observed correlation function 
£ (obs) (0), for which one obtains 



*(obs) 



d^rfe-^f'^i/)), 



(21) 



(22) 



where the self-convolved PSF T is given by 

T(9) = J d 2 iPW(Q-4>)W(lp); 

as can be easily seen, the normalization of W implies that T 
is also normalized. 

In analogy to the definition (8), one can now define the ten- 



sor M 



(obs) 



of second moments of the observed auto-correlation 



function; using Eq. (21), one finds that 



(obs) 



M -, 

where 



= Mij 



I 



d 2 e T{0)0i9j 



(23) 



(24) 



is the tensor of second moments of the self-convolved PSF. 
Hence, provided the PSF can be obtained with sufficient accu- 
racy, the tensor Mij can be directly obtained from (23), and 
thus the distortion 8 and magnification can be determined 
from (15) and (17) as before. The accuracy with which Mij 
can be obtained depends of course on the accuracy with which 
T(lp) is known, as well as on the relative 'size' of M and T; 
if the components of T are much larger than those of M, the 
determination of M = M^° bs:i — T involves subtraction of two 
'big numbers' and is thus subject to large fractional errors. 

If the PSF is isotropic, so is T; in this case, we can then 
write for the tensor %j = TSij. The distortion and magnifica- 
tion in terms of ( obs ^ is then given by 



Si 

82 

M 



trX(° bs ) - 2T ' 

2M ( ° 2 hB) 
trX(° bs ) -2T' 

VdetX(° bs ) - TtrX< obs ) + T 2 
M ' 



(25) 



2.3. Practical estimate of the distortion 
2.3.1. Estimates neglecting a PSF 

The calculation of the tensor M of second moments of the 
correlation function according to (8) is impossible in practice, 
since the integral extends over the whole 1R 2 , and therefore 
will be dominated by noise. One would like to have the inte- 
gral extend over a finite region only. Here, we outline several 
possibilities how, despite this difficulty, the distortion S can be 
determined. 

Adding a weight factor u>e(|0|) in the integrals in the defi- 
nition (8), 
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Jd 2 ege) e t e ]We (\e\) 

jd*8ao)M\o\) ' 



(26) 



would destroy the transformation relation (9), since this 
isotropic weight function would correspond to an anisotropic 
weight function in the source plane. However, the weight func- 
tion could be chosen to be about zero for small values 
of \6\ to minimize the PSF pollution, and also for large value 
when £ approaches the noise level. In particular, the distortion 
can be determined without any knowledge about the intrin- 
sic correlation function except that it is assumed to be 
isotropic. This is the method we choose for the practical im- 
plementation of the ACF method in the next Section. 

Another possibility is to define the tensor of second mo- 
ments as a function of the Jacobian matrix; for that we note 
that we can write the Jacobian as A = aB, where 



B 



1 + v 7 !- \5\ Z -5i -62 

-8 2 1 + V^ - W + SiJ ' 



(27) 



and a is a scalar. We now define the <5-dependent moment ma- 
trix 



Mi 



J d 2 6»£(6>) Orfj w (s) (\B0\ 
Jd 2 6>^(6>)w( s )(|B6»|) 



(28) 



where w^^x) is a weight function. Note that this matrix now 
explicitly depends on 8, and can be calculated from the ob- 
served correlation function for each value of 8. Changing the 
integration variable to (3 = AO — aBO, one finds that 



M {s) = AM A = a 2 BMB, 
where the matrix 

w /d 2 /3C (s) (m/W s) (1/3|A0 



M 



J d 2 /3£( s )(/3)wM(|/3|/a) 



(29) 



(30) 



is defined in analogy to (28). What is important to note is 
that the isotropy of ^ implies that M^ is again an isotropic 
matrix. In other words, the expression 



(M$-M%y+4(M t;) 

(a<<?+A<«) ! 



00 \2 



Y(sK2 



(31) 



should vanish. If M^ is replaced by the observed tensor M 
using Eq. (29), v will not vanish in general, unless the correct 
value for 8 was chosen. Performing this substitution, we find 
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«/(*) 



4detX (l - \S\ 2 ) 



[tiM - 5i(Mii-M 2 2) - 282M12] 



(32) 



= E 



[C (obs) (6> fc ) -Jd 2 t/> T(0 fc -^)£ W (|A^|)] : 



(33) 



note that the scalar factor a drops out. Thus, v depends on 8 
explicitly, as well as through the dependence of M on 8. Since 
v is positive semi-definite, a minimization of v will determine 
the value of 8. 

To summarize the basic idea of this approach, the weight 
function in the definition of the second moment tensor is chosen 
such as to yield an isotropic weight factor in the source plane, 
and it is adjusted such that the tensor of second moments in the 
source plane, as calculated from the observed tensor and the 
transformation (31), becomes isotropic. The weight function 
will be chosen such that it decreases to zero at radii where 
the signal-to-noise of the observed correlation function reaches 
about unity; for example, could be a Gaussian, with width 
adjusted to the data quality. 

2.3.2. Inclusion of a PSF 

If the effects of a PSF are included, together with a weighting 
scheme, the transformation becomes much more complicated 
than the one in (23). Hence, as was the case for the estimate of 
the distortion from galaxy images, the analysis becomes con- 
siderably more complicated if the effects of a PSF are taken 
into account (BM; Kaiser, Squires & Broadhurst 1995). The 
most serious effect is the gradual deterioration of the 'isophote' 
shape toward the ACF center. An issue to this problem is to 
use simulated images to calibrate the PSF effect on the shape 
of the objects as described in BM. We choose this method con- 
jointly to the biased weighting sheme (Eq. 26) instead of the 
unbiased Eq. (28). The reason is that for real data, the ACF is 
generally too noisy to allow a good shape parameter estimation 
from the intensity weighted formula given by Eq. (28) . 

Since the signal-to-noise of the observed correlation func- 
tion £( obs ) can be quite large, compared to single faint galaxy 
images, it might be worthwile to use an 'image restoration tech- 
nique' to invert the relation (21) between the observed corre- 
lation function £( obs ' and the true one, £. From the result of 
this restoration process (which could be, e.g., the Lucy 1994 
algorithm), the tensor M. can be calculated and the method 
outlined above can be used. Obviously, the result will not only 
depend on the quality of the data £( obs ), but critically on 
the knowledge of the (local) PSF, or more precisely, its self- 
convolved form T. Unfortunatly, the image restoration tech- 
nique has not given satisfactory result in the cases we have 
tested, and we will not pursue this idea further here. However, 
we think that it should be kept in mind and analysed in more 
detail. 

There might be a better way for the determination of the 
distortion (and the local magnification) in the presence of a 
PSF. The correlation ^ s - ) is due to very many faint galaxy im- 
ages if the brighter objects are removed (see below). Therefore, 
one might hope that this intrinsic correlation function is 
a universal function, being the same in all directions. Such a 
source correlation function will be obtainable from the Deep 
HST survey. If that were true, one would have a 'standard 
source' at high redshifts, and it is no surprise that that would 
help the analysis considered here tremendously. To wit, in this 
case one could compare the observed correlation function £< obs ) 
with the one calculated from the PSF and the intrinsic corre- 
lation function i.e., one can define 



where the sum extends over all values (pixels) on which the 
correlation function £( obs ' is determined, and cr^ is the standard 
deviation of the value of ^ obs ) at 9 k . By minimizing x 2 , the 
matrix A can be determined. As has become clear from the 
discussion at the end of Sect. 2.1, k and 7 cannot be determined 
observationally, but the quantities which can be determined are 
8 and \fj,\. Hence, for the use in (33), one should parametrize 
the matrix A as 



A = 



1 



^2| M | (l - \S\ 2 + V^W) 

1 + ^1-\S\ 2 -S 1 -ja 

-8 2 1 + ^1-^ + 8, 



(34) 



where the symbol '=' implies that the right-hand side is equal 
to A only in the non-critical regime, but that in the critical 
regime, the equivalent form of A in terms of \n\ and 8 is the 
same as in (34). Thus, if the intrinsic correlation function is 
known, one can determine |p| and 8 by inserting the form (34) 
for A into (33) and minimizing % 2 with respect to \fi\ and the 
two components of 8. 

3. Practical estimation in the weak-shear limit, and 
error analysis 

3.1. The practical method 

3.1.1. Practical estimate of the correlation function 

Up to now we have considered a lens mapping which can be de- 
scribed by its linearized form (1) over a sufficiently large solid 
angle in the sky. In practice, if one observes a cluster of galaxies, 
the linearization of the lens mapping will be valid only locally, 
on scales small compared to the typical angular scale on which 
the surface mass density of the cluster varies significantly. In 
those cases, one is interested in obtaining local estimates for 
the distortion and the surface mass density, i.e., the average 
in (5) should be extended over a limited area only. For our 
task, in order to avoid boundary effects of Fourier transforms 
in finite size fields, we choose to compute the ACF from the 
direct method which is more CPU time consuming than the 
Fourier analysis, but more accurate. In addition, the observa- 
tional frame contains bright objects, such as foreground stars, 
cluster galaxies, or bright field galaxies which would dominate 
the obtained correlation function if not removed of the aver- 
age in (5). Taking these difficulties into account, we propose a 
practical estimator for the determination of the local correla- 
tion function. 

Let Li denote the frame on which data are available. The 
bright objects can be identified and 'cut out'. This is done by 
using the FOCAS package to select the bright objects and re- 
place them by a shaded area. It is important to note here that 
this shaded area does not bias the estimate of the ACF, pro- 
vided it is not included in the ACF computation. The selection 
criteria must be such that the bright objects are completely 
covered by the shade mask. The resulting frame after cutting 
out these patches will have a 'Swiss cheese' shape and will be 
denoted by U' afterwards. On this field U' , one then defines the 
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mean intensity (in this section, we drop the superscript 'obs' 
for the observed intensity and correlation function) 



(35) 



keu' 



where the sum extends over all the N' pixels in W , and the 
intensity deviations i k — I k — I. From that, the correlation 
function of brightness fluctuations at position Oo is defined as 



mo ) = 



E^W-e/2)^ + e/2)mv>-g„|) 
E^(W-0 o |) 



(36) 



where the sums in numerator and denominator are taken over 
all possible values of ip such that if) ± 0/2 both lie in W , and 
w(6) is a weight function, typically chosen to be a Gaussian 
of width A9. £(0; Go) is often named Energy of the image U' . 
Hence, the expression (36) defines the correlation function as 
the weighted average over products of pairs of points with sep- 
aration vector 6 which can be selected within U' , and the con- 
tribution to the average is weighted by the function w which 
depends on the separation of the midpoint of the pair from the 
point Oq at which the correlation function is calculated. 



Fig. 1. A characteristic auto-correlation pattern at the faint galaxy 
angular scale. The ACF is computed from a 256 X 256 superpixcl 
of the I band of Q2345. Note in the central part of the ACF a 
crux pattern. This is a possible consequence of the loss of Charge 
Transfcrt Efficiency and/or the shift and add procedure. 

An example of a correlation function is shown in Figure 
1. Only the central part of the ACF (40 x 40 pixels) is repre- 
sented. The rest of the image contains noise plus interference 
terms which correspond to the overlapping of a galaxy by an- 
other galaxy. These cross-correlation terms are irrelevant for 
our purpose. 

3.1.2. Second moment measurement 

Eq. (21) shows that the observed ACF is a convolution of the 
true ACF with the square of the PSF, which is \[2 times larger 



than the PSF for a gaussian shape. No information is lost since 
the true ACF is also v2 times larger than a typical galaxy size. 
However, this effect must be taken into account for the compu- 
tation of the moments of the ACF, in particular for the choice 
of the filter. As pointed out by BM, the usual way to estimate 
the second moments of an object within a limiting isophote is 
not ideal. The convolution by the PSF changes the shapes of 
the isophotes which become nearly circular in the center of the 
object with little information left on the gravitational shear. 
The best way to estimate Mij seems to use an annular weight- 
ing function wg(8) centered on the ACF. Note that, contrary 
to the case of individual galaxies, the center of the ACF is, 
by definition, perfectly known. Since most of the gravitational 
shear signal is extending between 1 and 2 seeing disks around 
the center of a galaxy, a natural choice for wg is: 



w e (0) 



(s/») 2 



(37) 



where s is the angular size of the seeing disk. The shear 7 = 
|7|e llp thus inferred must be corrected. Indeed, the magnitude 
7 1 is underestimated due to the use of the weighting function 
wg , and to the PSF convolution. 

The problem of this underestimation was discussed earlier 
by several authors, as in BM, Fahlman et al. (1994), and Kaiser, 
Squires & Broadhurst (1995). BM concluded that this effect 
can be corrected using simulations from which they calibrate 
the measured shear intensity to the true one. They showed 
that the measured shear intensity is simply proportional to the 
true one (provided that the weak shear approximation applies) . 
Fahlman et al. showed that the calibration constant can be 
found from a streching of the images and convolving them with 
a known and large isotropic PSF. Kaiser et al. proposed to use 
HST images, stretch them by a known shear, reconvolve them 
with a PSF typical for ground-based observations, and then 
compare the resulting value of the measured 'shear' with the 
input value to obtain the correction factor between true and 
observed shear. 

As in BM, we choose here to use simulated images for the 
calibration. However, HST images provide an elegant way to 
calibrate the effect, as quoted in Section 2.3.2, but this requires 
a deep "empty field", and will be done in further work. 

3. 2. Errors coming from the noise level of the A CF 

This Section is devoted to the error measurement coming from 
the noisy aspect of the ACF. This is an important issue: con- 
trary to standard methods where the shear is inferred from 
the averaged shapes of many galaxies, the ACF method gives 
ultimately the shape of only one object, the ACF itself. It is 
therefore crucial to properly quantify the various sources of 
noise. 

In the first subsection, the signal-to-noise ratio of the ACF 
(S/N) is expressed in terms of the galaxy number N s , the 
field size, and the S/N of the individual galaxies, called ao. 
The second subsection quantifies the dependence of the shear 
measurement on the S/N of the ACF. 

3.2.1. Dependence of the S/N of the ACF on the image pa- 
rameters 

To simplify the discussion, assume that the image contains N g 
galaxies having the same size, the same profile, and the same 
surface brightness. We can then write the surface brightness as 
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1(0) W + EkiihiO-Ok), where 7 s ky is the sky 

brightness, i n (9) are the brightness fluctuations, and the third 
term is the contribution from the galaxies, situated at positions 
8k, with brightness profile I g . Subtracting the mean / from 

1(0), one obtains that i(0) = h(0-O k )-I % + i„(0) = 

Efcii Ig(9 — dk) + in(9), where I g is the mean brightness con- 
tributed by the galaxies. Since the galaxies are faint, and we 
are far from the confusion limit, practically all photons can be 
attributed to the sky brightness, which implies that the noise 
in(9) is Poisson noise almost entirely due to the sky brightness 
photons, i.e., of ky = (i\(9)) = 7 s k y , if intensities are measured 
in photon numbers. In particular, this means that the noise is 
uncorrelated with the signal. 

Inserting the foregoing expression for i(0) into (36), one 
finds that 



moo) = 



5>(|t/>-0o|) 



inW-6/2) i n (lp + B/2) 

K 

+ i u (i})-e/2)^i^ + e/2-e k ) 
k=i 

+ i n (V + 0/2) ^ 4(V> - 912 - 9 k ) 



+ j2U^-9/2-e 3 )j2h(^ + 9/2-9 k : 



(38) 



3 = 1 



One thus sees that £ contains contributions from galaxies only, 
from products of galaxy brightness and the noise, and terms 
proportional to the square of the noise. Considering 0/0, 
then the fact that the noise in different pixels is uncorrelated, 
and that the noise is uncorrelated with the signal contributed 
by the galaxies implies that the expectation value of £ is de- 
termined by the final term in (38). This can be split up into 
terms with j = k, and in those with j / k. The former thus 
contain the ACF of the galaxy brightness profiles and its result 
is proportional to N g . The latter contains cross-correlations of 
different galaxies, and is proportional to iVg 2 times the two- 
point angular correlation function of galaxy positions. 

More precisely, if we define the auto-correlation function of 
an individual galaxy as 



£ g (0) = / d 2 tf 4(0 -0/2) 4(0 + 0/2) 



(39) 



and neglect the possible angular correlation of galaxy positions 
[i.e., the terms with j =fc k in the last term in (38)], then one 
finds that the expectation value (£(0; 0o)) of the ACF becomes 



{ ,(0;0 o) )^ y^-;;; l) , g (0), 

Li, w (\V-9o |) 



(40) 



provided the weight function w(9) does not vary on scales com- 
parable to the size of individual galaxies. Correspondingly, the 
noise in the ACF becomes 



[iz-(m9 )}] 2 ) = *t ky 



E*«> a (W-0o|) 



(41) 



For the special case of a rectangular top-hat filter with N x x N y 
pixels, centered on 0o, the corresponding signal-to- noise ratio 
of the ACF becomes 



S/N = 



MO) 



(42) 



where now N s is the number of galaxies contained in the top- 
hat filter. We thus see that the signal-to-noise ratio is propor- 
tional to the number of galaxies within the filter. This should 
be contrasted to the case where the shear is determined from 
individual galaxy images; in that case, the signal-to-noise in- 
creases only as the square root of iVg . This fact shows that the 
ACF method will be superiour to the standard method if the 
number (density) of galaxies is sufficiently high. Furthermore, 
S/N depends on the filter scale, or in other words, depends on 
the number density n of the galaxies, S/N oc n y/ ' N x N y . 

The S/N also depends on the flux of the individual galax- 
ies. This is a consequence of the fact that, by definition, the 
objects in the ACF are weighted by the square of their flux 
[see Eq. (36)]. Therefore the S/N increases as the square of 
the flux of the individual galaxies; that is as a\. 

In this simple discussion where we neglected the profiles 
and the sizes of individual galaxies, we can therefore express 
S/N as follows: 



N 



, , 256 x 256 Ar 2 



(43) 



The constant K can be determined from the profile and size of 
individual galaxies; here, we obtain it from the hypothesis that 
all the galaxies have the same size, profile and flux. They are 
modeled by a gaussian profile of diameter 5 pixels. We define 
the S/N of an object (ACF or galaxy) as the mean amplitude 
of the object in an annulus of radius 3 pixels with a width of 1 
pixel centered on the object, divided by the sky photon noise. 
As explained in Section 3.1.2, the inner and outer radii of the 
annulus are justified by the fact that we avoid the central part 
of the object which is strongly polluted by the PSF. 

With these definitions, and for a square-shaped two- 
dimensional top-hat filter w, we find that the S/N becomes 



S_ 

N 



35 



256 x 256 / N 



N X N V V600 



)(t) : 



(44) 



This law was checked on a large number of numerical simula- 
tions. 

For large N g , the linear dependence of S/N with JV g is 
a strong gain of the ACF method compared to the standard 
method of coadding images method (BM) for which the signal- 
to-noise increases as -J N g . A comparison of the S/N for both 
the ACF and BM method is shown on Figure 2. In the case of 
deep observations (large N g ), the ACF will provide an easier 
and more accurate shear estimate from the shape matrix of the 
ACF than the coadded images of galaxies. 

3.2.2. Effects of noise on the shear estimate 

Because the gravitational shear is determined from the 
anisotropy of the ACF, the sources of noise on the ACF previ- 
ously discussed affect the estimation of the shear. In particular, 
in the case of a single object, as for the ACF, the noise degrades 
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Fig. 2. Value of the signal-to-noise ratio for the ACF method 
(straight lines) and for the method of coaddition of images (lines 
proportional to the square-root), as a function of the total number 
of galaxies on the field N g . The signal-to-noise of individual galaxies 
<tq corresponding to the different line styles are indicated. 



the isophotes and randomly polarize the object. Consequently, 
in order to detect a significant non-zero shear value, the S/N 
ratio of the ACF must be sufficiently large, or, in other words, 
for a given value of S/N, only shear values greater than some 
threshold 70 can be detected with sufficient statistical signifi- 
cance. 

Figure 3 is a result of a series of Monte-Carlo experiments 
where for each of them we measured the mean and the disper- 
sion of the shear of 100 ACFs. Each ACF has a gaussian shape 
of 7 pixels diameter. Its true ellipticity 7t rU c is known, and the 
orientation is zero (i.e., we choose 7t rU c to be real). A noise pat- 
tern with a given dispersion and mean zero is added on each 
ACF; this fixes the S/N of the ACF. In each experiment, the 
100 ACFs differ only in their noise pattern. We have performed 
80 experiments, with 10 values of S/N, and 8 values of 7t ruc . 
The ellipticities are measured using the spatial filter (37) with 
a diameter s = 6 pixels. We also define the dispersion: 



Fig. 3. The upper plot is the relation between 7 b s and 7truc- The 
proportionality constant is about a factor 4, and it is caused by the 
spatial filter (37) present in the second moment matrix (26). The 
lower plot is the 1-sigma detection level of the shear as a function 
of the S/N of the ACF. 



0"7 = Y ( [(7obsl - 7truc) 2 + (7obs2) 2 ] ) (45) 

with 7o bs = 7obsi + i7obs2 = |7obs| e 2l¥ \ At a given value of the 
S/N, the threshold 70 is defined as a 1-sigma detection of the 
ellipticity: 

70 = °7 • (46) 

Figure 3a shows the relation between 7 b s , the observed 
ellipticity value of the ACF and 7 trU o, its true value. The ap- 
parent discrepancy is an effect of the spatial filter, and will be 
corrected with the calibration of the shear (see section 4.1). 
The quantitative dependence of 70 on the S/N of the ACF is 
plotted in Figure 3b. 



3.3. Other sources of errors 

3.3.1. Errors from the observational characteristics of the 
galaxies 

The ACF method is ideal if the galaxies have the same flux, the 
same size, and the same profile. Unfortunately, this is not the 
case. A distribution in size, magnitude, and profile overweight 
some galaxies, and can change the statistical properties of the 
ACF. In principle, this could change the value of the shear. 
The detailed analysis of these effects, and how they could lead 
to new constraints on the faint galaxies properties, is a scope of 
the forthcoming Paperll. Nevertheless we can already propose 
a solution to the flux- weighting effect. 

If the galaxies are distributed in flux, they would (if they 
had equal profile and equal size) contribute to the ACF in 
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proportion to the square of their flux; hence, the brightest 
galaxies could dominate the resulting ACF. However, this could 
be partly avoided in the following way: define 7(^0; m > rrij) 
as the complex shear as determined from the ACF, which in 
turn has been determined by cutting out all objects on the 
frame brighter than rrij, as described in Sect. 3. 1.1. Selecting 
a set of flux thresholds rrij, and a set of weighting factors Vj, 
1 < j < M, one can form the average 



<7(0o)> 



E 



i=i 



which is an unbiaised estimate of the shear at for all choices 
of Vj. Note that by an appropriate choice of the weight fac- 
tors, the contribution by relatively bright galaxies can be sup- 
pressed. One can now try to choose the weight factors in such 
as way as to minimize the dispersion of the resulting estimate; 
this will be the topic of a future work. 

A size and a profile distribution are expected to destroy 
the perfectly concentric shapes of the isophotes in the ACF. 
However, with the typical smoothing length we used (a Gaus- 
sian function with 0.8-1.5 arcmin radius), the galaxy number is 
sufficiently high on a 3-hours exposure at CFHT to make this 
effect negligible and to induce an almost noise-free and regular 
isophote pattern on the ACF. 

3.3.2. Errors from the sources polarization distribution 

The intrinsic ellipticities of galaxies induce statistical fluctu- 
ations on their mean ellipticity, i.e., on the shear estimated 
from that mean. If a e is the dispersion of the intrinsic ellip- 
ticity distribution [assumed, as always, that this is isotropic, 
so that the mean of the (complex) ellipticity vanishes], then 
this intrinsic ellipticity distribution leads to a dispersion o" sta t 
of the resulting shear estimate of 



shape matrix of the ACF, this deformation can be approxi- 
mately corrected, as in Eq.(9) of Section 2.1: 



Mt 



J 4instA4obs^4i nst 



(49) 



As mentioned by BM, this approximation works provided that 
the instrumental deformation is constant in direction, and does 
not exceed 5-10 %. 



(47) 3.4- Choice of the local smoothing length 



Cstat 



(48) 



The ACF is computed in a given smoothing window of the 
image where the distortion is supposed to be approximately 
constant. The characteristic size of this area is called the local 
smoothing length. Its typical scale must be optimized accord- 
ing to the limiting shear value and the typical spatial resolu- 
tion we want to achieve. First, the number of galaxies 7V g must 
be large enough for the statistical fluctuations of the resulting 
shear not to exceed the desired accuracy of the shear mea- 
surement. That is, <7 s tat of Eq.(48) must be negligable when 
compared to |7obs|- Second, |7obs| must to be larger than 70 
(see Fig. 3b), which means that the S/N of the ACF must be 
sufficiently high. Since it is related by Eq. (44) to the number 
of galaxies, JV g is also constrained. The choice of the smooth- 
ing length is a compromise between these considerations. If it 
is too large, the small scales of the shear are averaged out. If it 
is too small, the noise and the statistical fluctuations induce a 
false value for the shear. Table 1 summarizes some constraints 
on JVg in the light of Sections 3.2.1 and 3.2.2. These values are 
only indicative since an exact model for the faint galaxies (size, 
profile, magnitude) is needed to make realistic predictions (see 
paper II). Such a realistic model will modify the law (44) and 
changes JV g . 



Table 1. Lower limits on 7V g required to measure a shear 70 
in a superpixel of size 256 x 256 with faint galaxies. 



Thus, if the shear is determined from the ellipticity of indi- 
vidual galaxy images, this provides a limit on the possible ac- 
curacy of this measurement. The dispersion of the intrinsic 
ellipticity is easily determined from the data, since it nearly 
equals the dispersion of the observed galaxy images. 

An intrinsic ellipticity distribution will also introduce a dis- 
persion on the shear as estimated from the ACF; however, in 
this case the corresponding dispersion can be less securely es- 
timated from the data itself. However, since the ACF method 
takes into account information from much fainter galaxies than 
can be detected and measured individually, we can take <7 B tat as 
an upper bound on the dispersion introduced by the intrinsic 
ellipticity distribution. 

3.3.3. Error from instrumental effects 

The PSF is never circular, and an anisotropic correction of 
the shear estimate must be applied. The principal sources of 
anisotropics are known (see BM). It is possible to measure the 
anisotropy from stellar profiles and to model it by an instru- 
mental deformation matrix Anst as the gravitational deforma- 
tion matrix introduced in Section 2.1, Eq.(2). Ai nst does not 
contain the PSF, but only the instrumental deformation. By 
applying the inverse instrumental deformation matrix on the 
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N n (ao) is the minimum number of galaxies with an individ- 
ual S/N of (To required in an area of 256 x 256 pixels (square- 
shaped top-hat window) to have a S/N of the ACF sufficiently 
large to measure a shear 70. JVstat is the minimum number of 
galaxies required such that 70 is three times larger than cr sta t , 
assuming an intrinsic ellipticity dispersion of cr e = 0.1 [see 
Eq.(48)]. 

Moreover, the choice of the smoothing length could depend 
on the position on the image since the shear varies spatially. 
Practically, an iterative procedure for the choice of this length 
is hampered by the very large computational time needed. An 
easier way is to find the mean "best" smoothing length from 
different tests. However, it should also be pointed out that for 
the reconstruction of the surface mass density of a cluster from 
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image distortions (e.g., Kaiser & Squires 1993; Fort & Mellier 
1994; Seitz & Schneider 1995, and references therein) a spa- 
tially homogeneous dispersion on the observed shear appears 
to be optimal, i.e., as long as the shear estimate is unbiased, 
these reconstructions do not require that a shear different from 
zero is measured in each smoothing window. 

4. Applications to data 

Though the ACF method seems theoretically attractive, its 
efficiency has to be demonstrated on real CCD observations. 
In this section, we first show that the ACF method works on 
realistic simulated CCD images where a shear pattern has been 
introduced. In a second step, we analyse the shear field on real 
data corresponding to the two fields already studied by Bonnet 
ct al. (1993, 1994). 

4-1. Simulated images 

We have used simulated images similar to those shown in 
BM and where photon noises, CCD defects, optical and atmo- 
spheric artifacts, or flat field residuals can be introduced. The 
ACF method is tested both on an image containing a lens, 
and on an image with no lens. Thus we can check that the 
method does not introduce spurious shear, and that the sim- 
ulations provide an accurate calibration of the shear as well. 
Figure 4 shows the shear pattern as it is directly measured on 
the simulated image containing a lens. The smoothing length 
is a top-hat filter of size 50" which corresponds to about 30 
galaxies detectable in each smoothing window. The elliptic- 
ity distribution is a flat distribution between 0.6 and 0.8. The 
resulting shear profile, azimuthally averaged with respect to 
the lens center, is displayed on Figure 5a. The errors bars are 
dominated by the statistical fluctuations of the galaxy elliptic- 
ities; they are determined from Eq.(48) using N s ~ 30. At the 
bottom, the 3 measured points before the calibration (49) are 
plotted (see subsection 3.2.2). To first order (weak shear ap- 
proximation), the calibration is a multiplicative constant and 
does not depend on the value of the shear and the lens model. 
However, it does depends strongly on the PSF and the spatial 
filter used to compute the second moments of the ACF. The 
calibration is found from simulated sheared images which re- 
produce real image conditions on telescopes. The constant is 
the ratio of the observed (uncorrected) shear to the true shear. 
In Figure 5b, the shear profile after calibration of the unlensed 
image is shown. One sees that it is compatible with zero shear. 

4-2. Real data 

4.2.1. The shear around Q2345 

The observation of a weak gravitational shear in the field of 
the double quasars Q2345+007A,B was made by Bonnet et 
al. (1993) from deep CCD images obtained at CFHT. This 
observational evidence of the gravitational origin of this dou- 
bly imaged quasar was confirmed later by Mellier et al. (1994), 
Fisher et al. (1994) and Pello et al. (1996) who found an excess 
of galaxies centered on the shear pattern, and also close to the 
quasars. They revealed a fairly complex lens, probably com- 
posed of a distant cluster of galaxies and some bright nearby 
galaxies. 

We analyse here the B and I images of Q2345 used in Bon- 
net et al. (1993), and we refer to this paper for the observa- 



Fig. 4. Simulation of a 4 hours exposure at CFHT in the B band on a 
3'5 X 3'5 field. The seeing is 0'.'7 with no tracking errors. Galaxies are 
lensed by an isothermal sphere (<r = lOOOkm/s), with a core radius 
of 4" located at 200" left from the field center. The lens redshift is 
0.17 and the mean redshift of the sources is 1. The segments show 
the local orientation of the shear. Their length is proportional to the 
shear intensity. 



tional details on the fields. We discard the I-band image for the 
shear analysis because the ACF is dominated by a crux pattern 
(see Figure 1) whose origin is unclear, but probably associated 
with residuals from the addition of many shift and add images 
which has been moved only along the pixel directions of the 
CCDs. Some residuals could be due also to bad charge effi- 
ciency. These defects are damaging for the shear analysis but 
they also reveal the high sensitivity of the ACF method (this 
suggests that the ACF could be also useful for testing CCDs) . 

Figure 6a shows the shear pattern found from a set of 
magnitude selected objects which clearly confirms the pattern 
found by Bonnet ct al. (1993). 89 objects are detected in the 
blue magnitude range [26,27]. They are surrounded by a cir- 
cle of size 30 pixels, whereas the rest of the image is put to 
zero. Since the galaxies are detected here using magnitude cri- 
teria, we adopt a completely different strategy than outlined 
before. The ACF is computed only from these selected objects 
since we detect galaxies here. The consequence is a consider- 
able increase of the S/N of the ACF, compared to the case 
that the rest of the image was not put to zero. This is an 
interesting alternative of the coadding image in BM, since in 
the present approach, the centroid of the individual objects 
need not to be determined. However, the ultra faint objects 
are lost with this strategy. The objects are spatially ordered 
into 8 groups of 10 — 15 members, and the ACF is computed 
in each of these groups. Since the signal is highly significant 
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Fig. 5. (a) 1-dimcnsional shear profile from the simulation of Fig. 
4. At the bottom the uncalibratcd measure points are drawn. The 
theoretically expected shear profile is plotted as the dashed line, (b) 
The calibrated shear as measured from the simulation without a 
lens. The measured values are compatible with zero shear 

(S/N ~ 40), and the shear found is high (between 10 and 15 
%), the measurement error is dominated by the intrinsic el- 
lipticities of the galaxies. The probability that each segment 
(independently of the others) comes from a random orienta- 
tion of the galaxies (the ellipticity dispersion is a e — 0.2) is 
around 10%. As in the Bonnet et al. paper, the shear is not 
corrected for tracking errors and PSF convolution since there 
are not enough stars to evaluate these effects (initially these 
observations were not done for measuring weak lensing!). 

Figure 6b shows the shear pattern computed from the ob- 
jects fainter than m_g = 26.5. The brighter galaxies are masked, 
as described in Sect. 3. 1.1, and the ACF is computed in a gaus- 
sian window of radius 0'.9. Since there are many more objects 
compared to Fig. 6a, the S/N of the ACF is also high (~ 40), 
and the shear around the expected excess of galaxies responsi- 



ble for the image splitting of the QSO is recovered. At least 30 
objects are detectable in each smoothing window, which gives 
a lower limit of the confidence level per segment of 99%, assum- 
ing the same ellipticity distribution for the ultra-faint galaxy 
population as for the ones with ms = [26,27]. Remarkably, 
another shear pattern is detected with the same confidence 
level at the bottom of the image. Its superimposition with the 
I-band CCD image (Figure 6c) clearly reaveals that the sec- 
ond shear pattern is centered on a visible excess of galaxies at 
the top-right of the I image. At first glance, the shape of this 
condensation is roughly elliptical and shows two density peaks. 
From the number of galaxies we see, it is probably a rather rich 
cluster of galaxies which has been undetected previously. 

4.2.2. The shear around C10024 

The central region of the cluster 00024+16 and an external 
region located at 6' from the cluster center were analysed by 
BMF, also from deep images obtained at CFHT. Again, we 
refer to this paper for details about the observations. 

The shear analysis done with the ACF of the central part 
of CL0024 is displayed in Figure 7a. The Gaussian smoothing 
length is 40", and the galaxies used are fainter than mj = 26.5. 
The shape does not show differences with the BMF analysis 
but the coherence of the signal around the cluster center is 
stronger. The inner pattern clearly shows an elliptical shape 
oriented along the second bisecting line, indicating an elon- 
gated shape for the cluster along this direction. Furthermore, 
the innermost central region shows a significant signal in the 
direction perpendicular to the second bisecting line. This sig- 
nal was not reported by BMF, probably because of the lower 
signal to noise ratio as compared with the ACF. Along with the 
general orientation reported above, the central shear resembles 
what is it expected for a bimodal mass distribution. 

The shear analysis of the large external region is more inter- 
esting. The ACF method was used with 3 differents magnitudes 
cuts, mv > 24, my > 25, and mv > 26. The corresponding 
shear maps are plotted in Figures 7b-d. For each segment, 
the smoothing length was 100". The general pattern found by 
BMF is also recovered here. The brightest map (mv > 24) 
is really similar, and the probable perturbation also seems to 
be detected, at the same position, though with a lower signifi- 
cance level. The intermediate magnitude map looks similar to 
the former one, but the perturbation is no more visible. On 
the other hand, a significant signal appears at the top right 
corner (the northern side of the cluster). It could be an effect 
of a second clump. This signal increases significantly on the 
faintest map, and there is no doubt at all that it is real. Once 
we go further away from the cluster center, the signal become 
noisier than on previous maps, but on the other hand, close 
to the cluster center the intensity of the shear is even higher 
than on the brightest and intermediate magnitude shear map. 
It suggests that we may be less polluted by cluster galaxies at 
such faint magnitudes. 

5. Conclusion 

This paper presents a new approach to probe weakly distorted 
galaxies down to the background noise limit. Using the auto- 
correlation of the image enables us to increase the signal-to- 
noise ratio of the shear considerably, mainly for two reasons: 
(1) There is no more need to detect individual objects and 
to define individual centroids and individual shape parameters 
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Fig. 6. Shear patterns found with the ACF method in the field of Q2345+007. (a) Shear pattern from 89 galaxies with ms=[26,27]. (b) 
Shear pattern from the galaxies with tub >26.5. (c) The shear pattern from panel (b) overlaid onto the I-band image 
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Fig. 7. Shear maps of C10024+1654 obtained with the ACF. (a) The central region of C10024+1654 with objects fainter than m B = 26.5. 
The segments give the orientation and the intensity of the shear. The coherence is remarkable. Note also the inner shear pattern which 
seems oriented along the second bisecting line. The two large vertical segments at the bottom center of the image are artifacts generated 
after the merging of two CCD images which overlapped at this position, (b) External part C10024 with objects fainter than my =24. The 
overall shape is similar as in BMF. Note the faint signal close to the perturbation reported by BMF which confirms that we probably 
detected a clump of matter toward this direction, (c) Same as (b) for galaxies with my >25. The main difference to the previous one is the 
detection of a clear signal at the top right, (d) Same as (b) for galaxies with my >26. The shear field at the top right is also significant 
and shows that another clump of matter must be located outside the field of view northward from the cluster center. 
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for these galaxies. (2) Even the faintest galaxies which cannot 
be detected by standard algorithms provide information to the 
ACF. Assuming there is a very dense population of distant 
background sources, they improve considerably the signal-to- 
noise ratio. A direct consequence of this correlation method is 
its ability to give upper bounds on these ultra-weak galaxies' 
number counts (see Paper II). 

We have demonstrated that the ACF can be analysed as 
a single object whose geometry is related to the ACF of back- 
ground sources modified by the gravitational shear. The sec- 
ond moments of the ACF computed over the whole field pro- 
vides the magnification map, provided the ACF of the sources 
is known. Since the intrinsic ACF is due to very many faint 
galaxies, one can hope that it is a universal function, i.e., the 
same in each direction on the sky, so that it can be obtained 
from very deep HST images. The relation between the observed 
shear and the true shear inferred from the ACF can be deter- 
mined (e.g., by simulations), and the intrumental/atmospheric 
degradations can be corrected in a standard way, provided the 
PSF is known, but with a strong benefit of the high signal-to- 
noise ratio of the ACF with respect to individual galaxies. In 
our approach we have implicitly assumed that this population 
is at high redshift, and does not correspond to low-redshift 
galaxies with low surface brightness. 

The applications on simulated and real data are very en- 
couraging. The ACF method easily recovers the shear profile 
found by BMF with another and independent method. How- 
ever, as shown in the I-band image of Q2345+007, the draw- 
back is the need of a very high quality of image aquisition. In 
particular, the shift and add method must be as random as 
possible, preferentialy with no shift in columns and rows, to 
avoid privileged directions. 

It is also possible to use the ACF method using only se- 
lected galaxies. The practical method is to select the objects 
with a given selection criteria (magnitude, or colour) and shade 
everything which is outside of a circle centered on each ob- 
ject. The ACF restricted to these objects can then be com- 
puted and analysed as was demonstrated here. This method 
was mentioned earlier by Kaiser (1992), though he did not go 
into further developments about its practical implementation. 

We are aware that the ACF method is sensible to a flux, 
size, and profile distribution of the galaxies. Although we pro- 
posed a solution for the flux problem, the size and profile prob- 
lem are more difficult and will be adressed later. 

About the future developments, two main points must be 
kept in mind. First, the deep HST empty fields will help in the 
determination of the intrinsic correlation function. This leads 
to the simultaneous determination of both the magnification 
and the shear using the x 2 minimisation (Eq. 33). However 
the question whether these empty fields are really empty is 
still open, though given that a one-dimensional function should 
be determined from observations (i.e., the radial profile of the 
isotropic intrinsic ACF) , this does seem relatively easily possi- 
ble. Second, the advantage of the ACF to avoid the reduction 
of data process (galaxy detection, first and second moment es- 
timation) may be increased if a reliable deconvolution of the 
PSF on the ACF is possible. However, the proper difficulties of 
the deconvolution algorithms, in particular the uncertainties in 
the choice of the number of steps for the deconvolution process 
is still a source of problems. 

The linear dependence of the S/N on the number of galax- 
ies is a strong argument in favour of the ACF method for the 



analysis of large fields where probably a small shear, but co- 
herent over large angular scales, exists (e.g., Blandford et al. 
1991, Kaiser 1992, Villumsen 1996, and references therein). 
The large number of ultra-faint galaxies is well adapted to the 
measurement of small shear. 

Another question which should be investigated in future 
work is how the shear determined from the ACF method 
is combined with shear measurements with the 'standard 
method', i.e., from obtaining the average ellipticity of observed 
galaxy images. The combination should be made in such a way 
as to minimize the dispersion of the resulting shear estimate. 

Acknowledgements We thank H. Bonnet, P.-Y. Longaretti, C. 
Seitz and S. Seitz for stimulating discussions, S. Roques and 
P. Marechal for enlighting discussions about the correlation 
function and M. Dantel-Fort, J. Bezecourt and J.-M. Miralles 
for their helps during the data reduction. This work was partly 
supported by the "Sonderforschungsbereich 375-95 firr Astro- 
Teilchenphysik" der Deutschen Forschungsgemeinschaft (PS), 
l'lnstitut National de Sciences de l'Univers (INSU) and the 
Croupe de Recherche Cosmologie. 

6. References 

Bartelmann, M., Narayan, R., 1995, ApJ, 451, 60. 
Blandford, R. D., Saust, A., Brainerd, T., Villumsen, J., 1991, 
MNRAS, 251, 600. 

Bonnet, H., Fort, B., Kneib, J.-P., Mellier, Y., Soucail, C, 

1993, A&A, 280, L7. 

Bonnet, H., Mellier, Y., Fort, B., 1994, ApJ 427, L83 (BMF). 
Bonnet, H., Mellier, Y., 1995, A&A, 303, 331 (BM). 
Brainerd, T., Blandford, R. D., Smail I., 1996, ApJ, in press. 
Broadhurst, T., Taylor, N., Peacock, J. A., 1995, ApJ, 438, 49. 
Dalcanton, J. J., 1996, SISSA preprint. 
Cole, S., Treyer, M.-A., Silk, J., 1992, ApJ 385, 9. 
Fahlmann, G. C, Kaiser, N., Squires, C, Woods, D., 1994, 
ApJ 437, 56. 

Fischer, P., Tyson, J. A., Bernstein, C, Guhathakurta, P., 

1994, ApJ, 431, L71. 

Fort, B., Mellier, Y., 1994, A&AR, 5, 239. 

Kaiser, N., 1992, ApJ, 388, 272. 

Kaiser, N., Squires, C, 1993, ApJ, 404, 441. 

Kaiser, N., Squires, C, Broadhurst, T., 1995, ApJ, 449, 460. 

Mellier, Y., Dantel-Fort, M., Fort, B., Bonnet, H., 1994, A&A, 

286, 701. 

Miralda-Escude, J., 1991, ApJ, 380, 1. 

Pello, R., Miralles, J.-M., Le Borgne, J.-F., Picat, J.-P., Sou- 
cail, Bruzual, C, 1996, submitted. 
Schneider, P., Rix, H. W., 1996, ApJ, submitted. 
Schneider, P., Seitz, C, 1995, A&A, 294, 411. 
Schneider, P., Ehlers, J., Falco, E. E., 1992, Gravitational 
Lenses, Springer. 

Seitz, C, Schneider, P., 1995, A&A, 297, 287. 
Seitz, C, Kneib, J.-P, Schneider, P., Seitz, S. 1996, A&A, sub- 
mitted. 

Smail, I., Hogg, D., Yan, L., Cohen, J. C, 1995, ApJ, 449, 
L105. 

Squires, G., Kaiser, N., Babul, A., Fahlmann, G. G., Woods, 
D., Neumann, D. M., Bohringer, H., 1996, in press. 
Van Waerbeke, L., Schneider, P., Mellier, Y., Fort, B., in prepa- 
ration (Paper II). 

Villumsen, J., 1996, MNRAS, in press. 



